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Abstract. A self-organized branching process is introduced to describe one-dimensional ricepile model 
with stochastic topplings. Although the branching processes are generally supposed to describe well high- 
dimensional systems, our modification grasps some of the peculiarities present in one dimension. We find 
analytically the crossover behavior from the trivial one-dimensional BTW behaviour to self-organized 
criticality characterised by power-law distribution of avalanches. The finite-size effects, which are crucial 
in the crossover, are calculated. 

PACS. 05.65.+b Self-organized systems - 05.70.Jk Critical point phenomena - 45.70.-n Granular systems 



1 Introduction 

Since thepioneering work of Bak, Tang and Wicsenfcld 
(BTW) [Q0, the sandpile model became one of proto- 
type abstract models exhibiting self-organized criticality 
(SOC). The original BTW model as well as its variants 
(see e.g. ||[I],[|,|f0]) consists of a cellular automaton slowly 
driven by stochastic perturbation. The state of each site 
is described by number of grains on top of it. (Actually, 
this number is rather the slope than the height, if we like 
to interpret the model as a real sandpile. However, in ID 
models, investigated here, the description through slope 
and height variables are strictly equivalent.) If the number 
of grains exceeds a threshold, the site becomes active, a 
toppling occurs and grains are transferred to neighbouring 
sites, which then may become active and the process con- 
tinues. The driving consists in adding grains on randomly 
chosen sites. The critical state is reached asymptotically 
in the limit of infinitely slow driving |8|] . Fully determin- 
istic versions were also studied, showing periodic |||l0| or 
self-similar but non-random behaviour [ pd| . 

Even though experiments on real sandpiles did not 
confirm SOC behaviour, due to inertia effects |12| , [l3| ,[l4], 
|l5l, |l6l , |l7i|lS|| , in the experiments using rice p9|,p0|| instead 
of sand it was found that large aspect ratio of the rice 
grains can lead to SOC behaviour p9[ , contrary to the 
case of sand, which has grains much closer to spherical. 

Another difference between a typical sandpile and ri- 
cepile experiments is that the ricepiles used in the ex- 
periments are quasi one-dimensional p9j,p0|. While the 
original BTW model in one dimension is trivial, there are 
several variants of ID BTW model which exhibit non- 
trivial behaviour §,0,E|,E|,|1|1E|]. Also the sandpiles 



on quasi one-dimensional stripes were investigated p6[ . 
Several one-dimensional models devised es pec ially for mod- 
elling the ricepiles were studied (^,|^,[|||^|3^||,||||||, 
p5j, pqp7|| . The models taking into account a possible long- 
range rolling of grains are able to describe the transition 
from SOC behaviour typical for ricepiles to the inertia- 
dominated behaviour of sand heaps 1 38 , 39]. 



Besides numerous exact results and renormalisation- 
group calculations (to cite only a few items of a vast bib- 
liography, see j|igl|l.il,i3llS)' the mean-field ap- 



proximation [}47f 48 



was very useful in clarifying the 
nature of the SOC state, even though it cannot give cor- 
rect values of the exponents below the upper critical di- 
mension. 

It was realised soon that the mean-field approximation 
for sandpiles is related to critical branching processes [|o| 
bjj . This idea lead to the introduction of self-organized 
branching processes ^,|||,|5J.|5^,^.^j, which describe 
the approach to the critical state. Similar approach consist 
in mapping the sandpile to the percolation on Bethe lattice 

& 

The approximation is based on the observation that 
in high dimension activity returns to the same site with 
very small probability. So, we can suppose that in each 
step the toppling occurs at a site, which never toppled be- 
fore during the same avalanche. Each toppling is mapped 
to one branching. Statistical properties of avalanches are 
determined by the probability p of branching. This proba- 
bility is itself determined self-consistently. If the avalanche 
is sub-critical, it does not fall off the system and average 
number of grains, and thus p, increases. If, on the other 
hand, the avalanche is super-critical, it surely falls off the 
system, which leads to decrease of the average number 
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of grains and decrease of p. It was shown [B2| , that this 
process sets the p exactly to the critical value, where the 
avalanche sizes s have power-law distribution P(s) ~ s~ r 
with mean- field exponent r = §. 

The purpose of this work is to modify the self-organized 
branching processes in order to describe one-dimensional 
ricepile models. Our model will be designed to comprise 
the one-dimensional BTW model as a special case. Clearly 
we cannot obtain correct values of the exponents. Our 
main question will be, whether there is a sharp transi- 
tion from trivial ID BTW behaviour to SOC behaviour 
or what is the nature of the crossover from the former to 
the latter. 

The paper is organised as follows. In the next section 
we define our version of the branching process, suitable for 
treating the one-dimensional ricepile. We find the condi- 
tion for the criticality and investigate the crossover from 
the trivial one-dimensional BTW behaviour to the critical 
branching process. The self-organization toward the criti- 
cal state is investigated in the section 3. We first define the 
self-organized branching process, then find the fixed point 
of the dynamics and show that it exactly corresponds to 
critical branching process. We finally investigate the in- 
fluence of finite size effects and find the finite-size scaling 
form. The section 4 concludes and summarises the work. 



2 Branching process for one-dimensional 
model 

2.1 Ricepile model 

The ricepile models were already thoroughly investigated 
by numerical simulations. In fact, there are two variants 
of the one-dimensional ricepile model. The so-called "Oslo 
model" |3y,[3l],|32,|33| supposes that the critical slope de- 
pends on space and tinre, and assumes new randonr value 
after each toppling event. Another approach |27],E8LE9| 
assumes that the toppling occurs with certain probabil- 
ity, which depends on actual slope. It is the second ap- 
proach, which we will follow in this article. It may be 
also noted that a two-dimensional model which also im- 
plements stochastic topplings was studied before ]59| . 

We recall shortly the definition of the model. We con- 
sider a chain of L sites. The state of site i, i = 1,2, ...,L is 
described by a slope z% — hi — /ij+i where the height hi is a 
non-negative integer, with boundary condition ft-L+r = 0. 
If the pile is in a stable state, a grain is dropped on the site 
i = 1. The update then proceeds for all sites in parallel. 
We look for all sites which satisfy at least one of the two 
conditions (i) it just toppled, (ii) its right-hand or left- 
hand neighbour toppled p7[ . If i is such a site, it topples 
with probability 1 if z% > 2, with probability a £ [0, 1] if 
z,i = 2 and with probability if Zi < 2. A toppling at the 
site i nreans that Zi is decreased by 2 and z^_i and Zj+i 
are increased by 1. 

For a = or a = 1 we recover the standard one- 
dimensional BTW sandpile model with critical slope z c — 
1 or z c — 2, respectively. In the intermediate region, < 



a < 1, self-organized criticality was found in numerical 
simulations, with avalanche exponent r = 1.55 ±0.02 (f29[ . 
However, it is not clear, what is the behaviour of the model 
for a close to either 1 or 0. It seems, that for a finite 
system the behaviour is SOC (modified by finite size ef- 
fects) only if a is not too close to 1 or p^,|60[. The 
behaviour of the system when the system size diverges 
and a stays close or 1 was not clarified. We would like 
to study this question within the approximation provided 
by a self-organized branching process. 



2.2 Characteristic functions 

From the technical point of view we will use the method 
of characteristic function (discrete Laplace transform), de- 
fined for a function f(s) on integer numbers s as /(£) = 

Er= cm. 

We will see that the distribution of avalanches have 
generic form 



P(s) 



-s/s 



(1) 



for large s. In the mean-field approximation or in the 
branching process we have r = 3/2, while in one-dimensional 
BTW sandpile the exponent is r = 0. The process is crit- 
ical, if the cutoff avalanche size sq diverges, sq — > oo. 

In the language of characteristic functions the behaviour 
(Q) translates in the properties of the singularity in P(£). 

Generally we have P(C) ~ (C — Co)'' + nonsingular part. 
For the one-dimensional BTW process we have r\ = — 1, 
while true branching process has r\ — 1/2. The cutoff is 
given by the distance of the singularity from the point 
C = 1, namely sq ~ l/|Co — 1|- The process is critical, if 

Co = l- 

We will also see that the characteristic function for the 
branching process is typically the solution of a quadratic 
equation. The singular part of the characteristic function 
comes from the square root of the discriminant -D(C) of 
the equation, i.e. P(C) ~ V^KC) + nonsingular part. 
Therefore, n — 1/2 and the cutoff is given by the solution 
of the equation Z?(Co) = 0. If D(l) = 0, we have s = oo 
and the process is critical. 



2.3 Branching process 

Let us first recall how the branching process is used to 
describe the simplest case of the sandpile model, for which 
in each toppling event two grains are transferred to two 
randomly chosen nearest neighbours (Manna model |B). 
There are Nq sites in state z — and N\ sites in state z = 
1. The branching process starts by dropping a grain to a 
randomly chosen site. The probability of becoming active 



(to topple) is p ■■ 



JVi 



Two new branches arise from an 



AT +JVi ' 

active site. Each of them is active with probability p and 
a tree is created iteratively. The branching process stops, 
when no active sites are present at the end-points of the 
tree. The number of active sites, or number of branchings, 
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corresponds to the size of the avalanche. The probability 
distribution of avalanche sizes can be easily obtained with 
the use of characteristic functions |£^,^3||5J,^,q6| and 
gives the mean-field value of the exponent r = 3/2 

Approximating the sand- or ricepile models by branch- 
ing process is well justified in high dimensions, where the 
activity returns to the same point with very small proba- 
bility. It seems, therefore, that the use of branching pro- 
cesses in the opposite limit, in one dimension, lacks sense, 
because the return of activity is very frequent. However, 
we can use a very simple property of the return of activ- 
ity to make the approximation sensible. Indeed, the most 
frequent case when the activity returns to the same site is 
described by the following process. 

If the site i is active (it topples) , a grain is transferred 
to site i + 1 which can become active. If that happens, 
another grain is transferred back to site i (and also to 
site i + 2, but it is not important now) and thus the site 
i may become active again. This observation leads to the 
modification of the branching process suitable for the one- 
dimensional case. We should take into account explicitly 
the return of the activity just in the next step. We will 
do it by setting different branching probabilities for site 
which was active just one step before (it is the site to the 
left) and for the site which did not have to be (the site to 
the right). 

Because the grains are added only on the site i = 1, we 
have Zi > Vi. The condition that the site topples with 
probability 1 if z > 1 ensures that Zi < 2 Vi. We denote iV a 
number of sites with z — a. So, picking randomly a site, 
we have probability p a — No/(Nq + N\ + N2) of having 
z = a, where a = 0, 1, 2. 

Let us now describe the construction of the branch- 
ing process corresponding to the one-dimensional ricepile. 
There are three types of the points on the tree created 
by the branching process, according to the value of z e 
{0,1,2}. We denote q a the probability that a point with 
z = a do branch. The points with z = do not branch, 
i. e. go = 0, while the points with z — 2 always branch, 
so qi = 1. The points with z = 1 branch with probability 
a, i. e. q\ = a. The approximation consists in supposing 
that if a site did not topple in the previous step, it has 
probability p a of having z = a, while if the site did topple 
in the last step, the probability of having z = a is modified 
due to the previous toppling to the value 



a) 



Pa 



q a+1 p a+1 

l^b=0 1b+l Pb+1 



(2) 



where we used p^ = q$ = for convenience. 

If a branching occurs at a site, two new branches ( "left" 
and "right") emanate from it. The probability that the 
right branch ends with a point with z = a is p a , while 
for the left branch the probability is p' a . This way the tree 
corresponding to the branching process is created. The 
above described rules are illustrated in the Fig. |lj. 

The root of the tree should be treated separately. The 
reason is that in the ricepile model the avalanche starts 
by dropping a grain always on the left edge of the pile, i. 
e. on the site i = 1. If it topples, it transfers a grain only 





=»- 




Pi 



1-a 



aft 



P 



b) 





Fig. 1. Illustration of the branching process. In a) the pro- 
cesses following a grain drop are depicted. Original configura- 
tions and their probabilities are in the left column, final ones 
in the right column. The possible final configurations resulting 
from a toppling are framed together with their non-normalized 
probabilities. In b) the correspondence is shown between one 
branching event and the toppling, in which one new grain is 
added and two grains (shaded) are displaced to the left and to 
the right from the toppling site. In c) a sample realization of 
the tree is sketched. The full circles placed on the right-hand 
branches correspond to probabilities p a , while empty circles on 
the left-hand branches have modified probabilities p' a . 



to the right, while the grain going left falls off the sys- 
tem. If we translated this feature to the description of our 
branching process, the root would consist either of a single 
non-branching point, or a point with a single branch (the 
right one) emerging from it. However, we are interested 
in the regime of long trees, where the different behaviour 
of the root from the rest of the tree is irrelevant. So, we 
assume that in the branching process also the root obeys 
the same rules as all other points. Thus, all points, includ- 
ing the root, have either zero or two branches emanating 
from it. 

The key quantity will be P^is), the probability that a 
tree consisting of n levels starting with a point of type a 
contains s branchings. The probability of having s branch- 
ings (i. e. avalanche of size s) is then P n (s) — J2 a PaP^(s)- 
We can easily derive the recurrence relation for P°(s) 
which becomes particularly simple if we use the charac- 
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teristic function. We obtain 



2 

E 

b,c=0 



4°(C) - (1 - Qa) + Qa( V P bP ' c PLl(OPn-l(0 • (3) 



A straightforward calculation leads to the following 
equations for the characteristic functions 



and 



Pn(0 = 1 

P r i(C) = l-a + aP n 2 (C) 



P n (s) = (api + p 2 )P%(s) for s > 



(4) 



(5) 



Therefore the basic quantity of interest will be the char- 
acteristic function P 2 (C) ■ All properties of the branching 
process can be computed from it. The set of equations 
(0) thus represent a single recurrence equation for P 2 (s), 
which in the limit n — » oo leads to quadratic equation for 
the stationary distribution P 2 (() = linin^oo P^(C)- We 
obtain explicitly 

^ 2 (C) = 

(a Pi + (1 - a)p 2 ) (1 - api - Pa) 
Pi + api 



ap 2 + 2 a (1 - a) p x p 2 + (1 - 2a) p 2 + Pi a 



2^,2 



P2 + api 



«p 2 (P 2 (C)) 2 



■P 2 (C) 



(6) 



length. Both one-dimensional BTW and branching pro- 
cesses have the same general form (pi) of the distribution 
of avalanches for large s, but the one-dimensional BTW 
process is characterised by the exponents r = 0, rj = — 1. 
Therefore, together with checking the criticality condition 
(0) we must also look at the behaviour close to the singu- 
larity. 



We will prove in the section 3.2 that in the thermody- 
namic limit our ricepile model self-organizes so that the 
parameters stabilise at values 



. 2a- 1. 

px =max(0, ) 

a 

p 2 =1 - a . 



(8) 



If we insert these values into the criticality condition 
(0), we find that it is satisfied for any value of a, including 
the limit values of and 1. At the same time we find that 
the singularity is always located at Co = 1 • (Indeed, as we 
discussed in section 2.2, the criticality of the process is 
equivalent to the condition Co = !•) However, we find that 
the type of the singularity corresponds to the exponents 
rj = 1/2, r = 3/2 (critical branching process) only for a's 
within the open interval (0, 1), while at the points and 1 
the model corresponds to one-dimensional BTW process. 
This can be easily interpreted in the language of sand- 
and ricepiles. Indeed, for a = and 1 the system recovers 
the behaviour of one-dimensional BTW sandpile, which 
does not exhibit critical behaviour in the usual sense. (In 
fact, the avalanche distribution does exhibit a power-law 
distribution: all avalanche sizes have the same probabil- 
ity, which corresponds to the power with exponent 0. But 
this is not the situation we usually describe as critical be- 
haviour.) 



2.4 Criticality 

The discriminant D(Q of the equation (0) depends on the 
parameters p\,p 2l and a. The branching process is critical 
if D(l) = 0. This implies the following relation 



-api - (1 -a)p 2 



2ap 1 p 2 +p 1 2 a 2 +p 2 2 = (7) 



which determines a surface in the parametric space. On 
this surface the process is critical and the distribution of 
avalanche sizes has a power-law tail with exponent r = 
3/2. 

However, the latter statement is not strictly true in 
the sense that if the coefficient at the quadratic term in 
the equation (|6|) is zero, the process is not a true branch- 
ing process, because each parent can have at most one 
offspring. It corresponds to a process with an exponential 
distribution of avalanche sizes, which we will call, in this 
work, a "one-dimensional BTW process" . The important 
feature which makes it different from a generic branching 
process is that there are no true branching points. Indeed, 
there may be a non-zero probability that the process stops 
at a given point, but there is zero probability to be split 
into more than one branch. Therefore, the process does not 
generate a tree-like structures, but linear chains of random 



2.5 Crossover behaviour 

The question arises, how the behaviour with exponent r = 
3/2 inside the interval [0, 1] crosses over to the exponent 
t = at the edges. As the critical behaviour is related 
to the singularities of the characteristic function, we will 
turn to the investigation of the function P 2 (C) in more 
detail. 

Indeed, we find that if we expand the solution of Eq. 
(ph for small values of the parameter p defined as 



p(0 = 



2a(l-a) 



(- 1 - l + 2a(l-a) 



(9) 



we can express the solution in terms of p and expand in 
the lowest order (for p 2 Cl) 



p2 (0 = --\- 2 -i 
p V p 



p(C) 



(10) 



While, as noted earlier, the exact solution for P 2 (C) has 
always the singularity of the type rj — 1/2 for C — > Co = 
1, the approximate behaviour (10) has a singularity with 
i] = —1 located at the point Co = (1 — 2 a (1 — a))~ J > 1. 
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When a goes to either or 1 , the value of Co approaches 1 . 
This suggests the following scenario. For large avalanches, 
i. e. 1 — C^Co - 1 the singularity at Co = 1 is relevant and 
the avalanche size distribution has a power-law tail with 
exponent r = 3/2. 

For shorter avalanches, i. e. 1 — C larger or comparable 
to Co — 1 the singularity at Co becomes dominant. There- 
fore, for short avalanches we have one-dimensional BTW 
behaviour P(s) ~ exp(— s/sq) with a cutoff 



«o = |i-Col' 



1 - 2 a (1 - q) 

2a(l-a) 



(11) 



The next step is investigation of the behaviour of sq 
when a approaches either or 1. We find it by expand- 
ing the expression for Co as a function of a around the 
points and 1, respectively. To make the notation more 
compact, let us introduce the variable /i G {0, 1}, which 
distinguishes the two limit points a — and 1. We can 
see from (fyl) that the cutoff diverges as 



1 



so 



\a — fi\ 



112) 



for a — ► n. 

On the other hand, sufficiently close to the singularity 
at C ~ ^ Co = 1 the exponent -q = 1/2 is relevant. The 
question is, how close to C = 1 one behaviour crosses over 
to the other. We have one-dimensional BTW behaviour 
for p 2 <C 1, while critical branching process behaviour for 
1 — p 2 <C 1. A typical crossover value Ccr can be found by 
solving the equation 



P(Ccr) = 2 



(13) 



The avalanche size distribution will exhibit the crossover 
around s cr = 1/|1 — Ccr|- F° r s "C s cr the one-dimensional 
BTW behaviour with exponential cutoff, diverging to in- 
finity for a — and 1, will apply, while for s 3> s cr the 
distribution will have power-law tail with usual mean-field 
exponent —3/2, and therefore exhibits self-organized crit- 
icality. 

The point of the transition between SOC and one- 
dimensional BTW when a approaches to 1 or lies in 
the diverging crossover value for the avalanche size. Simi- 
larly as in the case of s , by solving the equation (|13|) with 
defintition (0) we find the following limit behaviour 



1 



/■'■■I 



so 



(14) 



valid for s 3> 1 and a close to and 1 . The scaling function 
has the form F(x) ~ e~ x for i< 1 and F(x) ~ x~ 3 / 2 for 
x ^> 1. Indeed, we can find the Laplace transform of the 
scaling function as 



for a — ► /i. 

We can see, comparing equations ([l2| ) and (|l4|), that 
the cutoff for the one-dimensional BTW behaviour is asymp- 
totically equal to the crossover at which the critical branch- 
ing process behaviour sets on. This suggests the scaling 
form 



POO 



1 



F(- 



so(a) so (a) 



-) 



(15) 



-*(v-i) F (x) dx = y- vV-1 



(16) 



From here we obtain immediately the expression for the 
scaling function through the Bessel function of imaginary 
argument 



F(x) 



h(x) 



(17) 



The expected behaviour for x <C 1 and x 3> 1 can be 
verified directly by inspecting the asymptotic behaviour 
of the Bessel function. 



3 Self-organization 

3.1 Self-organized branching process 

In the basic setup of our branching process, all three pa- 
rameters a, pi, P2 are freely chosen. However, in the ri- 
cepile model the only free parameter is a. The number of 
sites with given z can change during an avalanche, so that 
also the probabilities p\ and pi are modified. This defines 
a flow in the space of parameters p\ , P2 ■ Our task now is 
to establish stable fixed points of this dynamics and check 
whether they satisfy the condition (Q). If that happens, 
we can conclude that the system is self-organized critical. 

There are four types of events, which can happen dur- 
ing an avalanche. Let us denote them T2, Tl, El, E0. In 
the event T2, the point with z — 2 receives a grain and 
topples. As a result, the number of sites with z = 2 is 
decreased by 1, JVjj — ► N<z — 1, and number of sites with 
z = 1 is increased by 1, N\ — ► N\ + 1. Similarly, in the 
event T\ point with z = 1 topples, Ni — > iVi — 1 and 
Nq — ► No + 1, in event El site with z = 1 receives a grain 
but does not topple, Ni — > JVi — 1 and A^ — ► N2 + 1, and 
finally in event SO site with z = receives a grain and 
does not topple, Ao — > Ao — 1 and Ai — > Ai + 1. 

Using the variables y G {T, _E} and a, 6 G {0,1,2}, 
let us denote s ay b,n the number of events of the type yb 
occurring at the level n within the branching process, 
on condition that the very first site had z = a. There 
are s ay b = S^o s ayb,n such events in the entire real- 
isation of the branching process. On average, there are 
(s-yb) = J2 a Pa( s ayb) events of the type yb. The averages 
(s y b) are of central importance for the dynamics of the 
self-organization and can be easily obtained as follows. 

For the characteristic function of the probability dis- 
tributionof the number of events s ay b, n we obtain an equa- 
tion analogical to (0). To study the self-organization, we 
will need only the average number of events, which is 
(s a yb,n), calculated as derivative of the characteristic func- 
tion. Hence 



Sayb.n) = qa^^iPc + p' c ) (s C yb, n-l) 



(18) 
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This is a set of three recurrence relations, which may be 
reduced to one equation only, by considering the relations 
(so y b,n) = and (si v b,n) = a {s 2y b,n}, valid for n > 1. If 
we take as the basic quantity the average (s2 V b,n), we get 
a recurrence relation deternrining a geometric sequence 



with quotient 



{S2yb,n+l) — K-(s2yb,n) 



ap2 + (P2 +api) 
p 2 + api 



(19) 



(20) 



We recognise in the stationarity condition k = 1 the equa- 
tion (J7J), implying the criticality of the branching process. 
Summation of the infinite geometric series gives imme- 
diately 



{Syb) = \Pb + (Pb+ p' b ) 



a Pi +P2 
1 



K 



(Sbyb,l) 



(21) 



where the initial conditions are given by (sbTb,i) — Qb and 
(sbEb,i) =l — qb- 

The self-organization of the branching process is due 
to the changes in the numbers N a , caused by the toppling 
(and non-toppling) events. These numbers determine the 
probabilities p a . Therefore, for fixed a the self-organized 
branching process (SOBP) S(a) consists of an (infinite) 
sequence of branching processes 



S(a) 



(22) 



.(°) J°) 



= {B(a,p^,p 2 u '),B(a,p^,p^),B(a,p^,p^),...} 

where B(a,pi,p 2 ) is the branching process determined by 
fixed parameters a,pi,p2, defined above. The branching 
processes within the sequence differ only by the values of 
the parameters p\ , and p2 ■ Let us consider the t-th branch- 
ing process in the sequence. When realised, it changes the 
original values of the numbers N a , or, equivalently, the val- 
ues of the parameters p a . The average change is uniquely 
determined by the average number of events (s y b) ■ So, 
the SOBP is entirely determined by the transition rela- 
tions connecting the values of the parameters in the t-th 
and t + 1-st step 



(23) 



(24) 



pf^-pf^T^lpf) 
for i e {1,2}. We find explicitly 

Ti(pi,p 2 ) = 

apx + (1 - a)p 2 - a(2 - a)p\ -p\- 2(1 - a)pip 2 

api + (1 - a)p2 + 2apip 2 +p\ + a 2 p\ 
T 2 (pi,p 2 ) = 

a(l - a)p\ + (1 - 2a)p\P2 
api + (1 - a)p 2 + 2apip 2 + p\ + a 2 p\ 



3.2 Fixed point 

The fixed point of the self-organization dynamics can be 
found immediately by equating the right-hand sides of 



equations (EJ) to zero. Direct solution of the two coupled 
equations gives three fixed points 



= o, 


P2=0 


(25) 


= o, 


p 2 = 1 — a 


(26) 


2a-l 
a ' 


p 2 = 1 - a . 


(27) 



The correct solution is determined by stability consider- 
ations. The relations (£4|) are linearised around the fixed 
points and the eigenvalues of the resulting matrices of rank 
2 are found. The result is that the fixed point (Eq) is al- 
ways unstable, while ( pq ) is stable for a £ [0, 1/2) and (27) 
is stable for a G (1/2, ljT For a = 1/2 the fixed points (26 ) 
and (p7| ) coincide and both of them are marginally stable 
(i. e. the eigenvalues have zero real part). 

Therefore, we find that the fixed point corresponds to 
the values of the probabilities 



i 2a-l, 

p\ = max (0, ) 

a 

p 2 =1 - a 



(28) 



which proves the already announced result of Eq.(H). 



3.3 Finite-size effects 

In the numerical simulations of the ricepile model p3| , |34| , 
f3(|,[30) attention is paid to the fact that the critical be- 
haviour is observed only for large enough systems with a 
not too close to neither nor 1. We have already shown 
how the crossover length blows up when a approaches the 
edge values or 1. It is obvious then, that for small sys- 
tems the crossover value of the avalanche size may not be 
accessible and the critical regime in the tail of the distribu- 
tion is not observed at all. In this subsection we will inves- 
tigate the consequences of finite length of the branching 
process. There are two phenomena where the finite size en- 
ters the problem. First, if the maximum number of gener- 
ations in the branching process is L, instead of infinity, the 
distribution of the avalanche sizes will not extend to infin- 
ity either, but will be bounded by s < s max = 2 L — 1. More- 
over, if we take for example pi = 1, p 2 = 0, a = 1, then all 
avalanches will have size L, therefore a peak at s = L will 
occur, P(s) — S(s — L). If we move slightly from this posi- 
tion by increasing p 2 and decreasing p\ and a, a structure 
of multiple peaks located at s = L, 2L — 1, 3L — 3, ... will 
appear. This makes the analysis very complicated. 

Second consequence is the shift in the self-organized 
value of the parameters p\ and P2, which for finite L will 
deviate from the critical values. Therefore, the avalanche- 
size distribution will develop an exponential cutoff in the 
form P(s) oc s~ 3 ' 2 exp(— s/si). 

As the first problem brings new particular difficul- 
ties, we will concentrate only on the second one. This 
makes the analysis less consistent, but feasible. Thus, we 
should stress that in the following we will suppose that 
the branching process in question has unbounded length, 
but the self-organization is made in such a way, that only 
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the first L generations of the branching process are taken 
into account. 

Instead of working with finite-!/ version of the equa- 
tions ( p3|) and (|24|), describing the approach to the fixed 
point, we can use the set of equations 



(sei) = (ST2) 
{SEo) = (sti) 



(29) 



which determine the position of the fixed point. The only 
information lost in Eq. ( p9| ) is the stability of the fixed 
points. However, we suppose the stability will not be af- 
fected by the finite-size effects. Therefore, we will rely on 
the stability analysis performed for L = oo also in the case 
of finite L and calculate the finite-size corrections starting 
withEq. @. 

The point is that the equations ( P9|) should hold also 
for finite L. In fact, the expression ( pl|) for the averages 
(s y b) assume the same form, only the factor (k— 1) arising 
in the L — oo version should be replaced by the factor 
K = (k — 1)/(k l_1 — 1). Assuming K small for large L, 
we can find p\ and P2 in lowest order in K . Then, we return 
to the definition of K and find that K oc L , confirming 
that our approach is consistent. 

Hence, for finite L we find, by solving the equations 
(p9|) to lowest order in l/L, for a G (0, 1/2) 



1 



Pi 



(2a -l) 2 L K L 2 > 



1-a ln(l- 7 ) ^ 1 



(30) 



and for a G (1/2, 1) 

2a - 1 5a 2 - 5a + 1 ln(l - 7) 



Pi 



V2 



a 



1 



(2a- l) 2 a 
1 - a ln(l - 7) 



2a -1 
where we denoted 

11 -2a 

7 



L 



L 



+ o{- 2 ) 



0{±) 



1 
•L 2 



(31) 



for a G (0, 1/2) 



2 1-a 

12a- 1 r 

7=- for a G (1/2,1) 

2 a 



(32) 



The above formulae confirm that the explicit limit L — > 00 
gives the same result as obtained previously when working 
directly with L = 00. 

Using these results we can find the position of the 
square-root singularity in the characteristic function for 
the avalanche size distribution, solving the equation 
-C(Co) = 0. The distance from 1 then determines the ex- 
ponential cutoff of the distribution. We find 



l/si = |Co-l| 



^+°(^) (33) 



where 



a(a) 



ln 2 (l 



■7) 



4a(l - a) 



for a G (0, 1) (34) 



and asymptotically for L — > cxd and a fixed the avalanche 
distribution becomes the function of s L~ 2 only, 



P(s; a, L) oc 7T 3 G(sL- 2 a (a)) 
and he scaling function has the form 



G(x) = x 



= -r-3/2 e -a 



(35) 



(36) 



This scaling holds well for all a with exception of the point 
a = 1/2, where we have 7 = and hence <r(a) = 0. Then, 
the next order in l/L takes over and the scaling changes. 
Let us use again the variable /1 G {0, 1}, which dis- 
tinguishes the two limit points a = and 1. The factor 
(j(a) diverges as u{a) ~ cto \ a — ^\ l f° r a - * A*, where 
(To = (ln2) 2 /4. Therefore, we can write the following scal- 
ing form for the avalanche size distribution 



P(s;a,L) <x L 3 \a-fi\ %G(sL 2 |a-/i 



co J 



(37) 



for a — » /i. 

We can see that the power-law distribution holds only 
for avalanches shorter than L 2 \a—/i\. In other words, if the 
parameter a is close to the end-points of the interval [0, 1], 
we need to have systems of the size L » i/\/l a — Ji\ for 
being able to observe any sign of self-organized criticality. 

In the above calculations we tacitly assumed that we 
beyond the reg ime we called "one-dimensional BTW" in 
the section 2.5. This means s 3> s cr . In fact, we can always 



reach this regime by choosing L large enough. Therefore 
the presence of the one-dimensional regime does not influ- 
ence the scaling behaviour for large L. More precisely, we 
should have L 2 \a—fi\ 3> s CI . But because s cr itself diverge 
for a — ► /i as Ja — /z| , we obtain a stronger condition for 
the scaling (B7p to be valid, namely 



L> \a- fi\ 



(38) 



if, 



//,. 



4 Conclusions 

We investigated analytically the self-organized critical ri- 
cepile model. We defined a self-organized branching pro- 
cess, suitable for one-dimensional problems. The model 
is characterised by the parameter a G [0,1], the proba- 
bility of toppling at a sub-threshold site. For both limit 
values a = and a = 1 the model is equivalent to the 
one-dimensional BTW model with trivial (uniform) dis- 
tribution of avalanches. 

We found that in the thermodynamic limit the sys- 
tem is self-organized critical for all values of a within the 
open interval (0, 1), with power-law tail in the distribution 
of avalanche sizes with mean- field value of the exponent, 
t = I. However, the power-law behaviour holds only for 
avalanches longer than a certain crossover value of the 
avalanche size. The crossover diverges when a approaches 
to either of the limit points of the interval [0, 1]. We found 
the scaling as well as the exact form of the scaling func- 
tion for avalanche distribution close to these limit points. 
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This describes how the one-dimensional BTW behaviour 
develops when approaching the limit points. 

The finite-size effects play important role in determin- 
ing whether the model is self-organized critical or not. 
In our model the SOC behaviour starts to occur at the 
larger sizes the closer we are to the limit points a = or 
1. We found the form of the finite size scaling in our self- 
organized branching process and determined the necessary 
condition for the the power-law regime in the avalanche 
distribution to be observable, when we approach to the 
limit points. 

I wish to thank Maria Markosova for numerous useful dis- 
cussions which motivated me to this work. I am indebted to 
Petr Chvosta for clarifying comments regarding stochastic pro- 
cesses. This work was partially supported by the Grant Agency 
of the Czech Republic, project No. 202/00/1187. 
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